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SUMMARY 


The  optimal  aerodynamic  design  of  new  configurations  requires  com- 
puter programs  which  will  calculate  surface  pressure  distributions  effi- 
ciently and  reasonably  accurately.  Numerical  methods  which  calculate 
the  entire  inviscid  flow  field  require  computational  times  and  storage 
too  large  to  use  in  design  optimization  processes.  Therefore,  a new 
method  has  been  developed  for  calculating  an  approximate  surface  pressure 
distribution  over  axisymmetric  bodies,  with  pointed  or  blunted  noses, 
at  angle  of  attack  in  supersonic  flow.  This  method  is  an  extension  and 
modification  to  the  second-order  shock-expansion  method  which  was  origi- 
nally developed  for  pointed  bodies  near  zero  angle  of  attack.  The  modi- 
fied Newtonian  pressure  distribution  is  used  in  the  nose  region  back  to 
a "matching"  point  where  the  local  Mach  number  is  slightly  supersonic. 

Then  the  new  method  is  applied  downstream  of  the  matching  point. 

In  the  late  60's  Jackson,  et  al , of  NASA  Langley  determined  a differ- 
ent matching  point,  but  it  did  not  pick  up  the  overexpansion  near  the 
juncture  at  sphere-cones  at  low  supersonic  Mach  numbers.  The  use  of  the 
new  matching  po'nt  predicts  this  overexpansion  quite  well.  One  modifica- 
tion made  to  thp  basic  second-order  shock-expansion  method  was  the  develop- 
ment of  an  "exact"  method  for  calculating  the  pressure  gradient  at  corners 
to  replace  the  approximate  method  normally  used.  Another  modification 
made  was  the  development  of  a new  method  to  calculate  the  effective  cone 
pressure  to  replace  the  tangent  cone  pressure  typically  used  for  bodies 
at  an  angle  of  attack.  This  effective  cone  pressure  was  found  to  be  much 
more  accurate  than  the  tangent  cone  pressure. 

The  new  method  has  been  applied  to  sphere-cones  and  sphere- tangent 
ogive-cone- flare  bodies.  Mach  numbers  ranged  from  1.5  to  4.63  and  angles 


xi 


of  attack  up  to  12°.  Surface  pressure  distributions  were  found  to  com- 
pare well  with  experimental  data.  In  addition,  the  resulting  axial  force, 
normal  force,  and  pitching  moment  coefficients  were  also  found  to  compare 
well  with  experimental  data.  The  computer  program  is  relatively  fast, 
requiring  about  15  seconds  for  a typical  case  on  the  IBM  370/165  computer. 
The  accuracy  and  simplicity  of  this  computer  program  makes  it  attractive 
for  engineering  applications. 
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INTRODUCTION 

Present  and  future  mission  requirements  of  the  Navy  dictate  a need 
for  methods  to  design  optimum  aerodynamic  configurations.  Of  particular 
interest  are  methods  which  will  predict  the  static  and  dynamic  aero- 
dynamic characteristics  of  wings,  bodies,  and  wing-body  combinations 
throughout  the  speed  regimes  of  subsonic,  transonic,  and  supersonic 
flow  and  at  angles  of  attack.  Since  a number  of  configurations  must  be 
analyzed  in  an  optimization  process,  it  Is  imperative  that  the  techniques 
employed  to  calculate  the  aerodynamic  coefficients  be  efficient  as  well 
as  accurate. 

There  are  generally  three  approaches  which  could  be  used  to  estimate 
aerodynamic  coefficients.  First,  wind  tunnel  and  ballistic  range  tests 
could  be  employed,  but  they  are  expensive,  time-consuming,  and  additional 
tests  are  required  for  new  configurational  innovations.  Secondly,  hand 
calculations  can  be  performed  using  handbook  techniques,  but  they  also 
require  a large  number  of  man  hours  and  their  accuracy  Is  questionable 
or  unknown.  The  third  method  is  to  develop  a computer  program  based  on 
a combination  of  analytical  techniques  and  semi -empirical  methods.  Al- 
though the  third  method  is  also  costly  and  time-consuming  initially,  It 
allows  the  designer  to  calculate  a large  number  of  configurations  and 
flight  conditions  accurately  and  expeditiously,  and  it  Is  ideal  for  de- 
sign optimization  techniques.  Therefore,  the  third  method  Is  considered 
the  best  approach  for  long-term  use  and  it  is  the  method  used  in  this 
report. 

A number  of  computer  programs  have  been  developed  to  calculate  the 
aerodynamic  characteristics  of  flight  vehicles,  but  none  of  them  covers 
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the  range  of  Mach  numbers  (0  < < 6),  angles  of  attack  (0  ± a <_  15°), 

and  complete  configurations  which  are  of  present  concern.  Quite  general 
body  and  wing  geometries  are  of  current  Interest,  and  they  include  spin 
stabilized  bodies  and  guided  and  unguided  fin  stabilized  bodies. 

Woodward1  used  a perturbation  method  to  compute  the  pressure  distribu- 
tion on  wing-body  combinations  at  subsonic  and  supersonic  speeds.  How- 
ever, his  method  does  not  include  blunted  bodies,  transonic  flow,  skin 

friction  drag,  base  drag,  and  nonlinear  angle  of  attack  effects.  Jackson 

2 3 

and  Smith  used  the  second-order  shock-expansion  theory  to  calculate  the 

drag  on  pointed  bodies  of  revolution  at  zero  angle  of  attack  in  supersonic 

flow.  They  also  calculated  the  pressure  and  skin-friction  drag  but  did 

4 

not  include  the  base  drag.  Jackson,  et  al  modified  the  second-order 
shock  expansion  method  for  blunt-nosed  bodies  and  angle  of  attack.  This 
method  was  found  to  predict  the  pressure  distribution  and  force  and 
moment  coefficients  on  blunted  cones  and  blunted-cones  with  flares  at 
Mach  numbers  from  about  1.5  to  5 and  angles  of  attack  to  12°.  Of  parti- 
cular importance  was  the  fact  that  their  adaptation  of  the  second-order 
shock  expansion  method  predicted  some  of  the  over  expansion  near  the 
sphere-cone  juncture  at  Mach  numbers  greater  than  about  2.0.  However, 
below  Mach  2 the  pressure  prediction  near  the  shoulder  was  poor. 

Another  method  for  calculating  aerodynamics  of  wing-body  combinations 

5 

is  that  due  to  Saffell  et  al  . It  is  only  applicable  to  low  aspect- 
ratio  configurations,  and  the  drag  calculations  are  inaccurate  at  small 
angles  of  attack.  An  empirical  method  called  the  "Spinner"  program^ 
was  developed  by  the  General  Electric  Company  for  spin-stabilized  bodies, 
however  it  will  not  handle  guided  vehicles.  There  are  several  computer 
programs,  e.g.,  Solomon  et  al^,  which  use  finite-difference  methods  to 
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calculate  the  entire  flow  field  around  a body.  However,  these  techniques 
are  generally  limited  to  inviscid  supersonic  flow  over  bodies  without 
wings  or  fins,  and  they  require  computational  times  much  too  large  to  use 
in  design  optimization  processes. 

The  best  computer  program  which  is  currently  available  to  compute 
static  aerodynamic  characteristics  of  general  wing-body  geometries  is 

g 

the  one  by  Moore  . This  program  predicts  the  pressure  distribution, 
lift,  drag,  and  center  of  pressure  for  Mach  numbers  from  zero  to  three 
and  angles  of  attack  to  fifteen  degrees.  The  results  compare  well  with 
experimental  data  and  other  analytical  results.  Also,  a typical  compu- 
tation costs  only  about  five  dollars  per  Mach  number  to  compute  the 
static  aerodynamic  characteristics  of  wing-body  configurations  on  the 

CDC  6700  computer.  For  the  pressure  distribution  over  the  body  alone, 

g 

Moore  employed  Van  Dyke's  second-order  perturbation  method  . Since 
this  method  is  designed  for  pointed  bodies,  Moore  used  the  modified 
Newtonian  pressure  distribution  near  the  nose  of  blunted  bodies  and 
then  matched  it  with  Van  Dyke's  second  order  perturbation  method  for 
the  region  downstream  of  the  nose.  This  procedure  was  shown  to  yield 
surface  pressures  which  compared  well  with  experimental  data  on  blunted 
cones  at  Mach  numbers  less  than  about  2.5  or  3.0.  In  particular,  it 
picked  up  the  over-expansion  of  the  pressure  near  the  shoulder,  whereas 
the  second-order  shock-expansion  method  was  found  to  be  more  accurate 
at  Mach  numbers  greater  than  about  2.5  to  3.0.  At  these  Mach  numbers 
there  is  no  over-expansion  near  the  shoulder. 

A two-year  contract  was  initiated  in  December  1976  at  the  North 
Carolina  State  University  to  extend  Moore's  method  to  Mach  numbers  from 
2.5  to  6 at  angles  of  attack  from  zero  to  15°.  This  report  describes 
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the  work  completed  during  the  first  year  of  this  contract  which  Is  re- 
stricted to  body^aTornr^er^dynaffilcs.  Modifications  are  made  to  the 
second-order  shock-expansion  method  to^frlculate  the  pressure  dlstrlbu- 
tlons  over  blunted  or  pointed  nose  bodies  with  flared  or  boattalled 
bodies.  The  skin-friction  and  base  drag  .are  computed  from  the  same 
techniques  used  by  Moore.  The  resulting  pressure  distribution,  skin 
friction,  and  base  drag  are  Integrated  to  yield  force  and  moment 
coefficients. 
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ANALYSIS 

Method  for  Predicting  Pressure  Distribution  at  o * 0 

Modified  Maslen  Method  - The  first  method  considered  for  calculating  the 
pressure  distribution  over  axlsymmetrlc  bodies  was  Maslen's  method^. 

This  method  was  modified  for  a = 0 and  used  by  the  first  author^  for 
blunt-nosed  bodies  In  hypersonic  flows.  Figure  1 compares  the  pressure 
distribution  computed  by  Maslen's  method  with  that  of  Solomon's  method7 
for  M^  * 3 and  o = 0.  This  figure  Illustrates  the  Inability  of  Maslen's 
method  to  pick  up  the  over-expansion  near  the  shoulder  at  the  lower 
supersonic  Mach  numbers.  At  hypersonic  speeds  the  over- expan si  on  does 
not  occur  and  Maslen's  method  predicts  the  pressure  distribution  very 
well.  It  Is  apparent  that  the  approximations  used  In  Maslen's  method 
are  not  accurate  at  the  lower  supersonic  Mach  numbers,  and  therefore 
this  technique  was  discarded  for  the  applications  considered  In  this 
project. 

Second-Order  Shock-Expansion  Method  - The  basic  second-order  shock- 

3 

expansion  method  was  developed  by  Syvertson  and  Dennis  for  pointed  bodies 
near  a = 0 at  supersonic  speeds.  In  order  to  apply  this  method,  the  ac- 
tual body  Is  replaced  by  a tangent  body  which  is  a series  of  conical 
frustrums  tangent  to  the  actual  body  at  a selected  number  of  positions. 
Figure  2 Illustrates  a typical  pointed  body  and  the  corresponding  tangent 

body.  The  pressure  on  the  Initial  cone  Is  obtained  from  a cone  solution, 

12 

e.g.,  the  method  of  Sims  . The  pressure  along  that  Initial  cone  Is 
constant,  but  the  pressure  drops  dlscontlnuously  across  the  juncture  of 
the  Initial  cone  and  the  following  conical  frustrum.  This  pressure  drop 
can  be  calculated  from  the  standard  Prandtl -Meyer  expansion.  In  the 
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generalized  shock-expansion  method  the  pressure  along  each  conical  frus- 
trum  is  taken  to  be  constant,  but  it  follows  an  exponential  variation  in 
the  second-order  shock-expansion  method.  This  exponential  pressure  varia- 
tion is  made  to  satisfy  three  boundary  conditions.  First,  the  pressure 
just  after  the  corner  of  the  initial  cone  and  conical  frustrum  is 
that  obtained  from  the  Prandtl -Meyer  expansion.  The  second  boundary  con- 
dition is  that  the  pressure  gradient  |||j  at  this  same  position  is  set 

equal  to  that  obtained  from  an  approximate  expression  developed  by 
Syvertson  and  Dennis  . For  the  third  boundary  condition,  the  pressure  at 
infinity  is  set  equal  to  the  cone  pressure  (pc)  that  would  exist  on  this 
conical  frustrum  if  it  were  infinitely  long.  Thus,  the  pressure  distribu- 
tion along  a conical  frustrum  is  given  by 


P * Pc  - (Pc  - P2)  e'n 


(1) 


where 


(x-x2) 

(Pc-P2)  cos  62 


(2) 


The  pressure  gradient  just  downstream  of  the  corner  is  determined  from 

. . 3 

the  approximate  expression 


where 


8 = -$4— 
2(M  -1) 


(3) 


(4) 


8 


'N 


X = 


2yp 

sin  2y 


(5) 


(6) 


Note  that  ||-|  is  the  curvature  of  the  surface  which  is  zero  on  conical 

frustrums,  and  fi  is  the  one-dimensional  area  ratio.  For  the  first  conical 

frustrum  after  the  initial  cone,  ||^-|  = 0 since  the  pressure  is  constant 

on  the  initial  cone.  The  pressure  gradient  ahead  of  subsequent  conical 
frustrums  is  obtained  from  the  derivative  of  Equation  (1). 

As  illustrated  in  Figure  2,  the  pressure  obtained  from  Equation  (1) 
is  used  only  at  the  position  where  the  conical  frustrum  is  tangent  to  the 
actual  body.  If  the  conical  frustrum  has  a negative  cone  angle.  Refer- 
ence 3 suggests  using  pc  = p^  on  that  segment.  Reference  3 also  gives 
an  expression  for  the  pressure  gradient  downstream  of  a sharp  concave 
corner  such  as  a flare.  One  restriction  on  the  method  described  here  is 


that  the  pressure  gradient  just  downstream  of  a corner. 


must  have 


the  same  sign  as  the  pressure  difference  ( PC_P2 ) • Otherwise,  Equation  (1) 
will  not  yield  the  third  boundary  condition  that  p = pc  as  x -*•  » . For 
more  details  of  this  technique,  see  Reference  3.  In  the  generalized  shock- 
expansion  method,  the  pressure  along  each  conical  frustrum  is  taken  to  be 
constant  and  equal  to  the  value  obtained  from  the  Prandtl-Meyer  expansion 
at  the  corner,  pg. 
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Jackson,  et  al  Modifications  to  Second-Order  Shock-Expansion  Method 

The  original  second-order  shock-expansion  method  was  developed  for  pointed 

4 

noses.  Jackson,  et  al  modified  it  for  blunt  nosed  bodies  by  using  the 

modified  Newtonian  pressure  distribution  up  to  a "matching-point"  and 

then  applying  the  second-order  shock-expansion  method  downstream  of  that 

point.  They  found  that  the  best  matching  point  corresponded  to  the  point 

on  the  blunted  nose  where  the  body  slope  is  the  same  as  the  maximum  wedge 

angle  for  an  attached  shockwave.  A tangent  body  is  used  to  replace  the 

actual  body  downstream  of  the  matching  point.  The  pressure  on  the  first 

conical  frustrum  tangent  to  the  matching  point  is  taken  to  be  constant 

and  equal  to  the  modified  Newtonian  pressure  for  that  body  angle.  The 

pressure  distribution  downstream  of  the  first  conical  frustrum  is  calcu- 

3 

lated  by  the  original  second-order  shock-expansion  method  . For  those 
conical  frustrums  which  have  the  sign  of  the  initial  pressure  gradient 

opposite  from  the  sign  of  (Pc-P2)»  the  pressure  is  taken  to  be  Pg. 

This  makes  the  technique  on  those  conical  frustrums  identical  to  the 
generalized  shock  expansion  method.  Reference  4 found  this  condition 
to  occur  on  the  nose  of  spherically  blunted  cones  at  > 3.0  and  also 
on  the  flares  of  blunted  cone-flare  bodies.  The  modifications  made  for 
bodies  at  an  angle  of  attack  are  discussed  below. 

Results  for  blunted  cones  and  blunted  cones  with  flares  were  found 
to  compare  well  with  experimental  data  for  ^3.0.  for  the  range  1.5  < 

M < 3 the  method  predicted  only  part  of  the  pressure  over-expansion  at 
the  shoulder.  However,  this  inaccuracy  was  found  to  have  only  a small 
effect  on  the  force  and  moment  coefficients. 


New  Modifications  to  the  Second-Order  Shock-Expansion  Method 

Two  modifications  are  made  here  to  the  modified  form  of  Jackson, 

4 

et  al.  . The  approximate  equation  for  the  pressure  gradient  given  by 
Equation  (3)  is  compared  to  an  "exact"  pressure  gradient  downstream  of  a 
corner,  and  a new  matching-point  is  determined  for  matching  the  second- 
order  shock-expansion  method  with  modified  Newtonian  theory  on  blunt- 
nosed  bodies. 

A.  Derivation  of  "Exact"  Pressure  Gradient  Downstream  of  a Comer 

As  mentioned  above,  the  pressure  gradient  given  by  Equation  (3)  is 
only  approximate.  Therefore,  an  "exact"  expression  is  derived  here  to 
compare  with  it.  Although  an  "exact"  expression  is  derived  in  Reference 
13,  it  is  not  in  a form  comparable  to  Equation  (3).  For  the  analysis 
below,  refer  to  Figure  3 for  the  nomenclature  and  geometry.  Consider  a 
surface  streamtube  as  it  passes  around  a comer.  The  stream  surface  (A) 
is  close  to  the  body  surface,  and  the  distance  along  a Mach  line  between 
the  two  stream  surfaces  is  "a"  (where  a is  very  small).  As  the  flow 
turns  around  the  corner  the  distance  "a"  will  increase  because  M increases 
and  hence  the  distance  between  the  stream  surfaces  will  increase. 

Define  C-j  and  C2  as  characteristics  coordinates  and  s and  n as  coor- 
dinates along  and  normal  to  a streamline.  At  a corner  only  the  family 
of  characteristics  (left-running  characteristics)  is  present.  The  analy- 
tical development  here  is  based  on  two  equations  from  the  method  of 
3 

characteristics  : 

. cos  » (f  - X f ) (7) 
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When  Equation  (7)  is  applied  at  the  corner  (on  the  surface),  there  results 
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where  the  subscript  "o"  refers  to  the  surface  streamline.  Since 


3C  i 

is  finite  and  the  distance  around  the  corner  is  zero,  dsQ  = 0 and  the 
right  side  of  the  equation  above  is  zero.  The  result  is  the  Prandtl- 
Meyer  expression 


dp 


d6. 
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and  since  dv  = - — 
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Equation  (9)  integrates  to 


v - v = 6 - 6 
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where  the  subscript  "o^1  refers  to  the  beginning  of  the  expansion  on  the 
surface.  Now  consider  the  streamline  off  the  surface,  A.  The  pressure 
and  deflection  angle  along  streamline  (A)  can  be  represented  by  a Taylor 
series  expansion  about  the  surface  streamline: 


Pa  = P, 
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The  partial  derivatives  are  evaluated  on  the  surface  from  Equations  (7) 
and  (8).  Since  the  distance  "a"  changes  with  6Q,  it  is  convenient  to 
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ratio  it  to  its  value  when  M = 1,  a*.  Then  it  can  be  related  to  ft , the 
ratio  of  the  one-dimensional  cross-sectional  area  of  the  streamtube  to 
that  at  M = 1 , as  follows: 

a _ Z-nra  _ ft  n.\ 
a*  ~ 2nra*  ~ sin  p ' ' 

Note  that  the  thickness  of  the  streamtube  is  a/sin  p and  thus  the  cross- 
sectional  area  of  the  streamtube  is  2Trra/sin  w.  Now  substitute  Equations 
(7),  (8),  and  (14)  into  (12)  and  (13)  to  obtain 
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Also,  use  the  Taylor  series  expansion  to  write 
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Appendix  A shows  that  the  differential  distance  along  streamline  (A)  is 
given  by 


(M2o-1) 


Substitute  Equations  (9)  and  (15)  through  (21)  into  Equation  (7)  applied 
along  streamline  A.  The  following  equation  then  results: 


d n , n sin  6 

ds.  g r 
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This  equation  can  be  integrated  numerically  around  the  corner  along  with 
the  Prandtl -Meyer  equation  to  determine  Qq,  and  using  Equation  (17)  the 


pressure  gradient 


|£j  is  obtained  from  Qq. 


It  should  be  noted  that  Equation  (22)  neglects  the  total  pressure 
3p. 

derivative  in  Equation  (19).  As  noted  in  Reference  3,  this  term 

is  generally  negligible.  To  be  certain  that  it  is,  the  accuracy  of 
neglecting  this  term  will  be  determined  during  the  next  year  of  this 
project. 

P 3 P 

Figure  4 compares  the  pressure  gradient  parameter  — , calcu- 

lated  from  the  approximate  method.  Equation  (3),  and  the  "exact"  method. 
Equation  (22).  The  pressure  gradient  is  presented  for  different  deflection 
angles  downstream  of  the  corner  of  an  initial  cone  with  a 9 half  angle 
for  M = 2.5  and  5.  This  figure  shows  that  errors  in  the  pressure 

OO 

gradient  parameter  as  high  as  35%  occur  for  M^  = 2.5,  and  this  maximum 
error  drops  to  about  27%  at  M^  = 5.  In  both  cases  the  "exact"  values 
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are  larger  than  the  approximate  values.  Figure  5 gives  as  a 

P« 

function  of  deflection  angle  downstream  of  a 40°  half-angle  cone.  At 

M =2.5  the  error  is  about  40%  at  S = 14°.  The  curve  for  M = 5 has 

00  00 

a different  feature  in  that  < 0 for  11°  < 6 < 40°,  and  errors 

Pco  3ti 

up  to  about  25%  occur.  Figures  4 and  5 illustrate  that  large  errors 
can  occur  when  using  the  approximate  relation  for  the  pressure  gradient. 
Note  that  since  the  initial  body  surfaces  are  cones  with  attached  shock 
waves,  3pt/3C^  = 0 and  Equation  (22)  is  the  appropriate  "exact"  relation 
for  these  figures. 

For  slender  bodies,  the  effect  of  using  the  exact  pressure  gradient 
rather  than  the  approximate  relation  is  small.  Figure  6 compares  the  re- 
sults from  both  techniques  with  the  pressure  distribution  from  Solomon, 
et  al7  for  a pointed  ogive  with  a boattail.  The  "exact"  pressure  gradient 
method  gave  slightly  better  results  for  the  boattail,  but  otherwise  both 
methods  were  very  close.  This  same  statement  cannot  be  made  for  thick 
bodies . 

B.  New  Matching  Point  for  Blunted  Bodies 

The  location  of  the  point  used  to  match  modified  Newtonian  theory  to 
the  second-order  shock-expansion  method  has  a large  effect  on  the  accuracy 
of  predicting  the  over-expansion  on  blunted  cones  at  the  lower  supersonic 

4 

Mach  numbers.  As  mentioned  above,  Jackson,  et  al  used  the  position 
where  the  body  slope  is  the  same  as  the  maximum  wedge  angle  for  an  attached 
shock  wave.  They  found  that  this  location  predicted  part  of  the  over- 
expansion at  the  shoulder.  The  results  compared  well  with  experimental 
data  at  2.30  < < 4.63,  where  the  over-expansion  is  small,  but  the 

comparison  was  poor  at  = 1.5,  where  the  over-expansion  is  large. 

Although  the  work  statement  in  the  present  contract  called  for  a method 
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applicable  to  the  2.5  < M < 6.0  range,  it  was  agreed  to  try  to  extend 

oo  ““ 

it  down  to  M = 1.2.  If  this  could  be  done,  then  the  subroutine  of 

00 

Moore's  computer  program  which  used  Van  Dyke’s  second  order  theory 
could  be  replaced  by  the  present  method  and  it  would  be  applicable  for 
the  full  supersonic  range  1.2  < Mm  < 6.0. 

The  first  modification  made  here  was  to  account  for  the  pressure 
gradient  on  the  first  conical  frustrum  at  the  matching  point  on  the  nose. 

4 

Jackson,  et  al  assumed  the  pressure  gradient  to  be  zero.  However,  both 
Equations  (3)  and  (22)  show  that  whenever  the  body  curvature  ( - 86/ as ) is 
discontinuous,  then  the  pressure  gradient  is  discontinuous,  even  if  the 
body  slope  is  continuous.  For  a blunted  nose  the  body  radius  of  curva- 
ture just  forward  of  the  matching  point  is  R-j , whereas  the  radius  of 
curvature  on  the  conical  frustrum  is  infinite  and  hence  ( 36/ as >2  * 0. 
Thus,  at  the  matching  point,  the  pressure  gradient  at  the  beginning  of 
the  first  conical  frustrum  is  given  by  both  Equations  (3)  and  (22)  as 


fEL-P  ** 
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(23) 


and  since  the  slope  is  continuous,  p2  = P-|  • The  properties  upstream  of 
the  matching  point  (subscript  "1")  come  from  the  modified  Newtonian 
pressure  distribution 

2 2 

p^  = pt  sin  6 + p cos  5 (24) 

Hence 
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where  R^  is  the  radius  of  curvature  at  the  matching  point. 
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The  pressure  distribution  downstream  of  the  matching  point  can  now  be 
calculated  by  the  second-order  shock-expansion  method.  Equation  (23)  gives 
the  pressure  gradient  on  the  first  conical  frustrum  (which  is  tangent  to 
the  nose  at  the  matching  point).  A tangent  body  is  then  transcribed  about 
the  remaining  body  and  the  technique  described  previously  can  be  applied. 
Note,  however,  that  the  pressure  gradient  at  each  corner  is  calculated 
from  the  "exact"  expression  given  by  Equation  (22). 

4 

It  was  found  that  the  matching  point  used  by  Jackson,  et  al  was  not 

the  best  one  when  the  pressure  gradient  is  used  on  the  first  conical 

frustrum.  Figures  7 and  8 show  the  pressure  distribution  over  a sphere 

computed  with  three  different  matching  points  (M  = 1.1,  1.2,  and  1.4) 

for  = 5 and  3.5.  The  results  are  compared  with  "exact"  inviscid 
14 

calculations  . The  pressure  distributions  with  matching  points  at 

M * 1.2  and  1.4  are  very  close.  Although  the  distribution  for  the 

matching  point  at  M * 1.1  differs  initially  from  the  other  two,  it 

blends  into  the  other  two  solutions  downstream  of  the  matching  point. 

For  the  lower  supersonic  range,  Figure  9 gives  the  distribution  of  the 

pressure  coefficient  over  a sphere  with  matching  points  at  M = 1.1,  1.15, 
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and  1.2  for  = 1.5.  The  results  are  compared  with  experimental  data  . 
Here  the  location  of  the  matching  point  is  more  critical  than  that  in 
Figures  7 and  8.  It  is  also  evident  that  the  modified  Newtonian  pressure 
coefficient  is  not  very  accurate  near  the  matching  points.  Although  an 
optimum  matching  point  cannot  be  determined  from  Figures  7 through  9, 
it  appears  that  it  lies  in  the  range  1.1  < M < 1.2.  Further  investigation 
is  needed  to  define  it  more  accurately  as  a function  of  M^.  Presently 

the  matching  point  Is  calculated  as  the  position  on  the  nose  where  the 
modified  Newtonian  pressure  distribution  gives  a local  Mach  number  of  1.15. 


Experimental  (Ref.  15) 


M = 1.1  ) 


M = 1.15 
M = 1.2 


Present  Method 
(M  is  Mach  No. 
at  Matching  Pt.) 


Modified  Newtonian 


Method  for  Predicting  Pressure  Distribution  at  a > 0 

An  accurate  calculation  of  surface  pressures  on  bodies  of  revolution 
at  a > 0 is  much  more  complicated  than  that  for  the  body  at  a = 0.  At 
a = 0 the  meridian  lines  on  the  surface  are  also  surface  streamlines. 
However,  at  incidence  the  streamlines  wrap  around  the  body  and  no  longer 
follow  the  meridian  lines.  If  the  shape  of  a streamline  can  be  deter- 
mined approximately,  then  an  equivalent  body  of  revolution  can  be 
generated  for  that  streamline  and  the  method  described  for  a = 0 can 
be  used  to  calculate  the  pressure  distribution  along  this  equivalent 


body  of  revolution  representing  the  streamline.  Several  streamlines 
must  be  used  in  order  to  get  the  circumferential  pressure  distribution 
along  the  actual  body. 

Jackson,  et  al4  assumed  that  the  surface  streamlines  in  the  windward 
and  leeward  planes  of  symmetry  could  be  approximated  by  the  body  shape 
obtained  by  transforming  the  body-axis  coordinate  into  wind-axis  coor- 
dinates. This  transformation  is  obtained  by  rotating  the  body-axis 
coordinate  system  in  the  plane  of  symmetry  clockwise  through  the  angle  a 
with  the  center  of  rotation  at  the  center  of  the  spherical  cap,  as  shown 
in  Figure  10.  The  windward  (e  = -90°)  and  leeward  (e  = 90°)  streamlines 
are  transformed  into  the  equivalent  bodies  of  revolution  represented  by 
the  body  coordinates  in  the  wind-axis  system.  (See  Figure  11).  On  the 
other  hand,  the  streamline  for  e = 0°  is  approximated  by  the  body-axis 
coordinates  of  a meridian  of  the  true  body  shape.  The  equivalent  body 
of  revolution  for  this  streamline  is,  therefore,  the  same  as  the  true 
body  at  a = 0.  Between  e = + 90°  and  e=0°,  the  streamlines  for  each 
radial  angle  e of  interest  are  assumed  to  follow  the  equivalent  bodies 
of  revolution  represented  by 
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( 


(26) 

(27) 

Figure  11  illustrates  these  equivalent  bodies  of  revolution  representing 
the  streamlines.  As  noted  in  Reference  4,  this  approximation  requires  a 
to  be  small  enough  for  the  stagnation  point  to  remain  on  the  spherical 
cap,  and  also  a must  be  small  enough  for  the  coordinates  of  the  equiva- 
lent body  in  the  leeward  plane  (e  = 90°)  to  remain  positive. 

For  the  analysis  here,  it  was  decided  to  modify  the  technique  des- 
cribed above  for  -90°  < 6 < 90°,  but  use  the  same  technique  for  the  wind- 
ward (e  = -90°)  and  leeward  (e  = 90°)  planes  of  symmetry.  One  reason  for 
the  modification  is  that  for  -90°  < 0 < 90°  the  stagnation  point  on  the 
equivalent  bodies  of  revolution  do  not  correspond  to  that  on  the  true 
body.  For  example,  at  0 = 0,  Equations  (26)  and  (27)  place  the  stagna- 
tion point  at  the  nose  (x  = 0,  y = 0)  of  the  actual  body  instead  of  its 
actual  location  for  a > 0.  Therefore,  the  method  used  herein  is  to 
assume  that  the  streamlines  on  the  spherical  cap  follow  the  meridian 
lines  of  the  sphere  in  the  wind-axis  system  up  to  the  sphere-cone  junc- 
ture. For  the  region  downstream  of  this  juncture,  the  streamlines  are 
assumed  to  follow  the  meridian  lines  in  the  body-axis  system.  Although 
the  approach  used  here  downstream  of  the  sphere-cone  juncture  is  also 
an  approximate  method,  the  streamlines  on  the  spherical  cap  are  the 
correct  ones  up  to  the  sonic  point. 

The  cone  pressure  used  in  Equation  (1)  by  Jackson,  et  al^  on  each 
conical  frustrum  was  the  tangent  cone  pressure.  For  an  axisymmetric 
body  at  a > 0,  the  tangent  cone  is  an  equivalent  cone  whose  slope  is 


Xw  = (Rn-x)(l-cosa)  /l-cos2a  + y si  not  si  no  + x 
Yw  = -y(l-cosa)  /l-cosza  + (Rp-x)  sina  sina  + y 
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the  same  as  the  slope  of  the  actual  conical  frustrum  with  respect  to  the 
wind  axis  system.  The  equivalent  cone  angle  6 corresponding  to  a cone 
angle  6 is  given  by 

sin  = cos  a sin  6 - sin  a cos  6 sin  e (28) 

It  was  found,  however,  that  even  for  a pointed  cone  at  a > 0,  the  tangent 
cone  pressure  differed  from  the  actual  pressure,  particularly  at  the 
lower  supersonic  Mach  numbers.  In  order  to  get  a more  accurate  pressure 
on  the  tangent  cone  at  a >0,  the  tangent  cone  pressure  was  replaced  by 
the  following  relation: 
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The  second  term  on  the  left  is  the  cone  pressure  at  a = 0 whereas  the  two 
terms  on  the  right  are  the  modified  Newtonian  pressures  at  a > 0 and  a = 0, 
respectively.  This  equation  gives  the  change  in  cone  pressure  with  a as 
the  change  in  modified  Newtonian  pressure  with  a.  Since  the  modified 
Newtonian  pressure  is  given  by 
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IpJ  Ip*  J e(i 


and  6 = 6 at  a = 0,  then  the  pressure  on  the  tangent  cone  at  a > 0 

becomes 
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To  test  the  accuracy  of  Equation  (29),  comparisons  were  made  with 

14  15 

"exact"  numerical  calculations  and  experimental  data  . Figure  12 
shows  the  circumferential  pressure  variation  on  a pointed  5°  cone  at 
M =3.5  for  a = 0,  5°,  and  10°.  At  « = 10°  the  tangent  cone  method 

oo 

overpredicts  the  pressure  about  15%,  whereas  Equation  (29)  is  less  than 
1%  above  the  "exact"  calculation.  Figure  13  compares  the  circumferential 
pressure  coefficient  determined  from  Equation  (29)  to  experimental  data 
for  a 9°  cone  at  a = 12°  with  M ' = 1 .9.  In  the  windward  plane  the 

CD 

tangent  cone  method  overpredicts  p/p  by  20%  while  Equation  (29)  under- 

oo 

predicts  it  by  3%.  Figure  14  compares  the  present  theory  and  tangent 
cone  method  with  "exact"  numerical  results  in  the  windward  plane  for  a 
9°  cone  at  M =3.5  for  0 < a < 10°.  Here,  also.  Equation  (29)  was  found 

oo  — — 

to  give  better  results  than  the  tangent  cone  method.  As  the  Mach  number 
increases,  the  difference  between  the  results  from  the  two  methods  de- 
crease. This  fact  is  illustrated  in  Figure  15  for  a 20°  cone  at  M =10 

oo 

for  a = 1°,  5°,  and  10°. 
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Figure  15.  Circumferential  Pressure  Distribution  Around 
Pointed  20°  Cone,  Mm=10,  a=T0°,  5°,  0° 


RESULTS  AND  DISCUSSION 


The  pressure  distribution  and  force  and  moment  coefficients  have 
been  calculated  and  compared  with  experimental  data  for  blunted  cones 
and  blunted  cones  with  a flared  afterbody.  Two  of  these  bodies,  called 
Model  1 and  Model  2,  are  illustrated  in  Figure  16. 

Blunted  Cones 

Pressure  distributions  over  Model  1,  a spherically  blunted  11.5° 
cone,  at  a = 0 are  shown  in  Figures  17,  18,  19,  and  20  for  M =1.5,  1.9, 

oo 

2.3,  and  4.63,  respectively.  The  results  from  the  present  method  are 

4 4 

compared  with  that  of  Jackson,  et  al  and  experimental  data  . Figure  17 
and  18  indicate  the  capability  of  the  present  method  to  calculate  the 
pressure  over-expansion  near  the  shoulder  much  better  than  Jackson's 
method.  As  M increases  the  over-expansion  decreases,  and  at  M = 4.63 

oo 

it  vanishes. 

Figure  21  compares  the  calculated  pressure  distribution  with  experi- 
15 

mental  data  on  a spherically  blunted  9°  cone  at  a = 0 and  M^  = 1.5. 

Two  numerical  solutions  are  presented,  one  with  a matching  point  at 
M = 1.10  and  the  other  at  M = 1.15.  The  over-expansion  here  is  signifi- 
cant, and  the  results  for  M = 1.10  are  very  close  to  the  experimental 
data.  Figures  22  and  23  give  results  for  the  same  body  at  Mm  = 1.9  and 
4.63,  respectively.  The  agreement  with  experiment  is  good  here  also, 
and  the  over-expansion  is  not  present  in  Figure  23.  The  "hook"  in  the 
calculated  pressure  distribution  at  the  matching  point  is  the  same  as 
that  observed  in  Figures  7,  8,  and  9. 


Figure  17.  Pressure  Distribution  on  Model 


Pressure  Distribution 


Pressure  Distribution 


Experiment  (Ref.  15) 


Blunted  Cone-Flare 


Pressure  distributions  over  Model  2 are  compared  with  Jackson's 
4 4 

theory  and  experimental  data  in  Figures  24,  25,  26,  and  27  for  a = 0 
and  M^  = 2.3,  2.96,  3.95,  and  4.63,  respectively.  Since  the  over- 
expansion is  small  for  these  Mach  numbers,  the  calculated  results  are 
close  to  Jackson's  results.  Both  methods  compare  well  with  experimental 
data  except  near  the  beginning  of  the  flare.  It  appears  that  boundary- 
layer  separation  occurs  upstream  of  the  flare  and  thus  the  theory  does 
not  calculate  the  correct  pressure  in  this  region. 

Effect  of  Angle  of  Attack 

Figures  28,  29,  and  30  show  the  pressure  distribution  on  the 
blunted  9°  cone  at  a = 12°  and  M =1.5  for  0 = -90°,  0°,  and  90°.  This 
example  is  a severe  test  for  the  theory  because  it  is  a small  cone  angle 
at  an  angle  of  attack  greater  than  the  cone  angle  and  a low  supersonic 
Mach  number.  The  results  given  in  Figure  28  are  for  the  windward  plane 
(0  = -90°).  Note  the  improvement  in  the  calculated  pressures  when 
Equation  (29)  (labeled  "new  soln  for  pc"  in  the  figure)  is  used  to  calcu- 
late the  equivalent  cone  pressure  (pc)  in  comparison  with  the  tangent 
cone  method  (labeled  "old  soln"  in  the  figure).  The  results  for  0=0° 
(Figure  29)  and  0 = 90°  (Figure  30)  (leeward  plane)  do  not  compare  well 
with  experimental  data  although  the  present  method  is  somewhat  more 
accurate  than  Jackson's  method. 

The  inaccuracies  in  the  surface  pressure  distribution  appear  to  have 
only  a small  effect  on  the  force  and  moment  coefficients.  Figure  31 
gives  the  coefficients  of  pitching  moment  (C^),  normal  (C^),  and  axial 
force  (C^)  as  a function  of  angle  of  attack  for  Models  1 and  2.  This 
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figure,  the  same  results  as  Jackson's  method,  and  the  modified  Newtonian 
theory  predicts  accurate  results.  However,  this  fact  cannot,  in  general, 
be  stated  for  other  bodies  over  the  Mach  number  range  of  interest. 


i 
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Pressure  Distribution 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  Three  new  modifications  have  been  made  to  the  second-order  shock- 
expansion  method  which  makes  it  more  accurate  in  predicting  the  surface 
pressure  distribution  over  pointed  and  blunted  bodies  of  revolution. 
These  modifications  are: 

A.  An  "exact"  method  for  calculating  the  pressure  gradient 
around  a corner. 

B.  A new  procedure  for  matching  second-order  shock-expansion 
theory  with  modified  Newtonian  theory  for  blunt  nosed  bodies. 

C.  A new  method  for  calculating  the  effective  cone  pressure 
which  replaces  the  tangent-cone  method  for  bodies  at  angle 
of  attack. 

2.  This  new  method  gives  results  which  compare  well  with  experimental 
data  and  predicts  the  over-expansion  near  the  shoulder  of  blunted  cones 
at  Mach  numbers  less  than  2.0. 

3.  The  new  method  has  been  found  to  yield  accurate  forces  and  moments 
for  pointed  and  blunted  bodies  at  1 .5  _<  M^  £ 4.63  and  0 <_  a < 12°. 
Pressure  distributions  are  also  accurate  near  the  windward  plane,  but 
tend  to  be  inaccurate  near  the  leeward  plane  of  cones  with  cone  angles 
less  than  the  angle  of  attack. 

4.  It  is  recommended  that  additional  results  be  calculated  and  compared 
with  experimental  data  to  determine  the  accuracy  of  the  method  at  Mach 
numbers  as  low  as  1.2.  Also,  modifications  to  the  angle  of  attack  solu- 
tion are  needed  if  accurate  pressure  near  the  leeward  plane  are  desired. 
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GLOSSARY  OF  TERMS 


f 


a 

B 

C, 


C-j  ,^2 
t, 

M 

n 


P 

pt 

Q 

r 

R 


s 

x 

y 


a 

Y 
6 
n 
e 
x 
u 

V 


distance  from  body  to  streamline  A 
function  defined  by  Equation  (4) 
axial -force  coefficient 
pitching-moment  coefficient 
normal -force  coefficient 
pressure  coefficient,  Cp  = 2(p-poo)/(pooVw) 
characteristic  coordinate 
body  length  (calibers) 

Mach  number 

coordinate  normal  to  a streamline  (calibers) 

o 

static  pressure  (dyne/cm  ) 

total  pressure  at  stagnation  point(dyne/cm  ) 

function  defined  by  Equation  (17) 

radius  of  body  of  revolution  (calibers) 

body  radius  of  curvature  (calibers) 

coordinate  along  a streamline  (calibers) 

coordinate  along  body  axis  (calibers) 

coordinate  normal  to  x in  plane  of  sytmetry, 
see  Figure  10  (calibers) 

angle  of  attack  (degrees) 

ratio  of  specific  heats  (1.4  for  air) 

flow  inclination  angle  and  cone  angle  (degrees) 

function  defined  by  Equation  (2) 

meridian  angle,  see  Figure  10  (degrees) 

function  defined  by  Equation  (5) 

Mach  angle,  u * sin"^(l/M) 

Prandtl-Meyer  expansion  angle  (degrees) 
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$ angle  on  sphere  measured  from  nose  (degrees),  see  Figure  7 

n ratio  of  cross-sectional  area  of  streamtube  to  that  at 

M=l,  see  Equation  (6) 

Subscripts 

1 condition  at  start  of  Prandtl-Meyer  Expansion 

2 condition  at  end  of  Prandtl-Meyer  Expansion 

o surface  streamline 

o-j  beginning  of  the  expansion  on  the  surface 

A streamline  off  the  surface 

c quantities  evaluated  for  cone  tangent  to  the  body 

eq  refers  to  equivalent  cone  quantities 

MN  modified  Newtonian  value 

N nose 

w refers  to  wind  axes 

» free  stream  conditions 

Superscripts 

* evaluated  where  M=1 
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